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1 Introduction 



Nowadays, that the Gribov ambiguities play an important role in the quantization of Yang- Mills 
theories is becoming more and more evident. These ambiguities, affecting the Faddeev-Popov 
quantization formula, deeply modify the infrared behavior of the theory and might have a crucial 
role in the understanding of color confinement. 

The aim of these notes is to give a simple introduction to Gribov's work and to its conse- 
quences on the infrared regime of nonabelian Euclidean gauge theories. The discussion will 
be focused mainly on the original paper by Gribov P, whose content will be reproduced in a 
somewhat detailed way. 

This work is divided in two parts. Part I is devoted to the study of the Gribov ambigui- 
ties. We begin by reviewing the examples provided by Gribov pQ. Further, a class of Gribov 
copies proposed by Henyey [2] will be considered. 

In Part II we introduce the so called Gribov horizons and we analyze their properties. Fol- 
lowing Gribov's suggestion, the restriction of the domain of integration in the Feynman path 
integral to the first horizon will be discussed. The ensuing modifications on the infrared behavior 
of the gluon and ghost propagators in the Landau gauge will be worked out. 

We conclude this short introduction by mentioning a result due to Singer 0. Although the 
examples of the Gribov copies which we shall consider in the following will be worked out in the 
transverse gauge, dA = 0, the existence of the Gribov ambiguities is not related to a specific 
gauge condition. As pointed out in [3], the presence of Gribov copies is in fact a general feature 
of nonabelian gauge theories. 



2 Part I: The Gribov Pendulum 
2.1 The Gribov ambiguities 

2.1.1 Quantization of Euclidean Yang-Mills theories. Non-uniqueness of the gauge 
condition 

The Euclidean Yang-Mills action* 

Sym = 4^2 / ^xF^ , (1) 

is left invariant if one replaces by the gauge transformed field A^ 

A il ^A il = S^d^S + SU„S . (2) 

The fields A^ and A^ define equivalent configurations, leading to the same value for the expres- 
sion In order to properly quantize the theory, only inequivalent configurations should be 
taken into account in the Feynman path integral. The first step towards the implementation 
of this procedure is provided by the Faddeev-Popov quantization formula. One integrates over 
configurations which have a certain divergence, i.e. 

= f. (3) 

*See Appendix A for the notation. 
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Therefore, for the partition function one gets 

Z = N J DA IM 5(dA - /)det (M ab (A^j e'^ 1 d ^ xF ^ F ^ , (4) 

where 

M ab (A) = -d^5 ab - f**A%) , (5) 

is the Faddeev-Popov operator and N a normalization factor. In the following we shall limit 
ourselves to the choice / = 0, which amounts to impose the following transversality condition 

dA = . (6) 

This condition is known as the Landau gauge condition. Thus 

Z=N J DA^5(dA)det (^M ab ^j e'^ 1 ^^ . (7) 

However, as pointed out by Gribov pQ, the condition d A = does not fix uniquely the gauge 
configurations. This means that, for a given A^ satisfying the gauge condition dA = 0, there 
exist equivalent fields A^ obeying the same condition, i.e. 

A„ = Sfld lt S + &A li S = A lt + rf(d li S+[A li ,S\), 

dA = dA = 0. (8) 

Field configurations A^ satisfying the conditions (jSJ) are copies of A„. In terms of the gauge 
transformation S, the condition dA = dA = reads 

d ll (S\d tl S+[A tl ,S}))=0, (9) 

which, due to dA = 0, becomes 

d^d^s + s^d^s + d^A^s+s^A^s = o . (10) 

• Remark 

At the infinitesimal level, S = 1 + a, a -C 1, expression (jiflj) reduces to 

d 2 a-{d ll a)A ll + A fl {d ll a) = 0, (11) 

i.e. 

d IM {d^a + [A^a}) = 0. (12) 

We see that, in the infinitesimal case, the condition for the existence of Gribov 's copies 
is equivalent to state that the operator —d^ (d^ ■ + [A^, ■]), whose determinant enters the 
Faddeev-Popov quantization formula Q ; has zero eigenvalues. Notice that the eigenvalues 
equation for the Faddeev-Popov operator —du {dn ■ + [A^, •]), 

-d ll {d^ + [A^]) = e{A)^ , (13) 

can be seen as a kind of Schrddinger equation, with A^ playing the role of the poten- 
tial. Therefore, for large enough values of the field A^, we might expect that zero energy 
solutions, e = 0, of eq. (|iffi do in fact exist. 
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2.1.2 The Gribov pendulum 



In order to analyze the condition for the existence of Gribov's copies, we shall consider the 
three-dimensional case, Aj, i = 1,2,3, assuming that the gauge group is SU(2) and that the 
gauge field is spherically symmetric, i.e. Aj depends on the unit vector n« = Xi/r, r = y/xiXi, 
riiUi = 1. Let <7j denote the 2x2 Pauli matrices 

o-jOj = Sij + ie ijk a k . (14) 



It follows that the quantity n 
obeys 



n = iriiOi , (16) 



n 2 = —ninj(Ti(jj = —UiUi = —1 . (17) 
Since the gauge field A, is Lie algebra valued, it can be parametrized as 

M = ai(r)Irii + a 2 (r)(Ji + a 3 (r)n i ("n > • ~a) + a±{r)ai(n ■ ~a) , (18) 



where I stands for the 2x2 unit matrix^ and where we have taken into account that A% is 
spherically symmetric. Notice that, due to rurii = 1, and ( n ■ ~a) 2 = 1, higher powers of (njnj) n 
and (n • ~~a) m are absent in the expression (|18|). Also, we observe that 

ieij k nja k = rij {oiOj - Sij) = <7i("n • 7r ) — n* , (19) 

meaning that this term is not independent. In addition, from eq. (|14j) . it follows that 

(n • ~o)(Ji = — o"j( n • cr ) + 2rij . (20) 

Thus, expression (|18|) yields the most general form for a SU(2) field which is spherically sym- 
metric. Moreover, due to the traceless condition 

Tra a A1 = , 

we get the relationship 

a 4 = -ai , (21) 
so that Ai turns out to be parametrized by three independent quantities 

Ai = ai(r)Irii + a 2 (r)ai + a 3 (r)n i ("n > • ~a) - ai(r)oi(~n ■ ~a) . (22) 

Of course, we can adopt the same parametrization employed by Gribov pP, namely 

fin f}fl 

Ai = fi(r)-£- + f2(r)n— + / 3 (r)nn, . (23) 



Indeed, from 



7— = -a k (5 ik - mn k ) = ~(<Ti-(n ■ a )th) , (24) 
OXi r r 



it follows 

i i _> _> 1 _> _> 1 
Ai = ~h{r)ai - ~/i(r)( n ■ <r )n» / 2 (r)( n ■ a )a { + -i^r)^ + i/ 3 (r)( n • cr )n» (25) 

y» 

which has precisely the same form of expression (1221) . 
^The set (I, en) is a basis for the 2x2 matrices of SU(2). 
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Observe that for fi = fs = 0, expression (p^) becomes 
A: = 



1 1 

-f2\r)(n ■ o-)oi + -f 2 {r)rii 



—h(r) + itjiko-k) rij + -/2(r)n, 



which is purely transverse, i.e. 



jSijkXj0- k f2(r) , 



<9,;A; = . 



Having found the most general parametrization for the SU(2) gauge field, ea. (|23|) . let 
out the condition for the existence of copies, i.e. 



where we shall consider the class of gauge transformations S parametrized by 



air) _> _> a(r) 
cos — h i n • a sin — — - 



From ens.fpHjl. (HSU it follows 



A; 



a 



cos ■ 



n sm • 



a ^ a\ 1 da a dn 

sin h n cos — - — h sin — — — 

2 2 / 2dxj 2 dxj 



(OL Oi\ I dfl 0TI \ / OL OL 

cos --resin -J {fi(r)—+f 2 (r)n— + f 3 (r)nnij (^cos- + resin - 



dn 



n 



1 9a 



a 



9re 



a^ dn 



re cos h re sin — - — h cos — sm — — sm" — re— — 

2 2 J 2 dxi 2 2 dxi 2 dx { 

( a dn a^dn a ^ a^dn 

+ I/ 1 C ° S 2 ^ + h C ° S 2 + /3 C ° S 2 nH * ~ h Sm 2 n dx~ x 

a dn a \ / a 

+72 sm- — + / 3 sm-rei ) ( cos - + /; sin - 



Since 

we obtain 
I; 



da , 
_ = a (r)ni 



1 ,^ 1 . <9re . 2 a ^dn 
— a (r ren, H — sin a— sm —re 

2 v ; 2 dxi 

, ( 2 a dn a . a 

+ji cos — — h cos — sm 



2 dx; 



a 



dn 



a 



2 ctej 
9re ^ 
dx^ 
a^ dn 



n 



2 a ^ dn ^ 
sm —re— — n 
2 dxi 

a a dn 



a 



dn 



+/2 I cos — re— h cos — sin — n— — re + cos — sin— — h sin — — — n 



2 dx 



2 dx 



2 dxi 



2 dxj 



+f$nni 



Observing that 



dn k 
dxi 



rij(a k o-j + o-jGk) 



„dn k 

2——n j d kj 

dxi 



dn k d(n k n k ) 
2^— n k = — 

dxi dxi 



0, 



one has 



and 



Therefore 



dx~ n = ~ n dx~- ' (34) 



^ dn ^ dn dn 

n^—n = -nn—— = —— . 35) 

oxi oxi axi 



I I 2 a 2 a \ /■ a a \ 3™ 

A = ( /i ( cos - - sm - J + 2/ 2 cos - sin - J — (36) 



, ,. , o Q ■ 2 a \ nJ - a . a\ ^dn 
+ 1/2 cos — — sm — — 2/1 cos — sm — n- 



2 2) J 2 2) dx 



Finally, recalling that 



I \ I dn S dn 

+ ( fa H — a (r) rim H — sin a— (1 — cos a) n— — 

1 J6 2 K ' 2 dxi 2 y ' dxi 



cos — — sm — = cos a , (37) 

a a 
2 cos — sin — = sin a , 
2 2 



we get 



Ai = I fi cos a + I f 2 + ^ J sin a 



2 / / <9x 



+ ( (/ 2 + 2~ C ° S °) ~~ ^ Sin ° ~ 2~) + (/ 3 + 5°'^) ™ Hi ' ( - 38 ' ) 

It remains now to work out the condition d{Ai = diA{. After a little algebra, one obtains the 
differential equation 



a "i r ) + ~ a 'i r ) ~ ^ {j^ 2 + ^ Sina + ^ ^ COSa " ^ = ' 



(39) 



Setting 

r = log r , r = e r , 

equation (|3*9*|) becomes 

a2a(T) +^-4ff/2 + ^sina + / 1 (cosa-l ) )=0, (40) 



<9t 2 <9t 

which is the equation of a pendulum in the presence of a damping term a'(r), see Fig.l. 



7 




Figure 1: The Gribov pendulum 

Summary 

We have started with the most general spherically symmetric 577(2) field 

Ai = fi( r )-Q^; + -MOn^ + h(r)nni . (41) 
For the gauge transformed field 

A % = S^dtS + S^AiS, (42) 



we have 



4 = ZiW^ + ^W^ + ZsW™.. (43) 



7i(r) = fx cos a + f / 2 + ^ sin a , 
h{r) = h + ^ c °s a - fi sin a - - , 
fs(r) = h + \a'(r). 

Finally, the condition diAi = d{Ai yields 

^1 + ^1 _ 4 ((/ 2 + sina + /i (eosa - 1)) = . (44) 

T/ifs is i/ie equation of a pendulum in the presence of a damping term a'(r). The compo- 
nents fi, $2 of the gauge field correspond to the forces acting on the pendulum, Fig.l, 
see also Appendix B. 
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2.2 Examples of Gribov's copies 

In this section we shall work out explicit examples of Gribov's copies. We shall restrict ourselves 
to the transversality condition 

d t A t = d t A t = . (45) 
As initial gauge configuration we shall take 

i 

M = -^eij k Xja k f(r) , diM = , (46) 

corresponding to setting f\ = = 0, fa = f in expression (|4T1) . The gauge transformed field, 
ea.(|42|). is given by 

f + 2 ) Sm ^ + [ f + 2 C ° S ~ 2 J "fl£ + 2° (r)nni • (47) 
Also, the condition = 0, gives 

a"(r) + -c/(r) x ( f + - J sin a = , 

r r z \ 2/ 

or, r = logr, 



(48) 



dr 2 8t V 2 
which corresponds to the damped pendulum of Fig. 2. 



d 2 a(r) da{r) A ( t 1\ . 

y + — ^-4 / + - sma = 0, (49) 




Figure 2: 
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The presence of the function f(r) in eq. lflfij) is needed to ensure that expression is a regular 
configuration. More precisely, we shall require that Ai is regular at the origin, r = 0, and that it 
goes to zero at infinity, r = oo. Essentially, we shall consider two types of boundary conditions, 
namely the so called weak (WBC) and strong (SBC) boundary conditions 0J|S]. 

• For the WBC we have 

Ai is regular at the origin, /(r) r _*o — ► 0(r) , 

Ai decays as 1/r for r — > oo, /(r) r _ >00 — > constant . 



(50) 



For the SBC 



Ai is regular at the origin, /(r) r _>o — ► O(r) , 

Ai decays faster than 1/r for r — » oo, /(r) r _ >oc — ► , 

(r^)™ ->0. (51) 

Let us begin with the case of SBC. Recalling that dn/dxi ~ 1/r as r — > oo, it follows that the 
equivalent field Ai in ea. ((17j) will obey SBC if 

a(r) — > 27rm + 7r /or r — > , m integer (52) 
a(r) — ► 27rm + 7e T /or r — > — oo , 

and 

a(r) — > 2-7rn /or r — ► oo , n integer (53) 
a(r) — ► 2irn for r — > oo . 

Two situations are possible, according to the strength and the orientation of the force /. 
• i) The first case corresponds to (1 + 2/ ) > 0, see Fig.3. 




Figure 3: 
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At r — ► — oo, the pendulum starts at a position of unstable equilibrium, a = 27rm, with 
velocity a' ~ 7e T , 7 / 0, see Fig. 4. 



1 1— ^ ce=2nm 



Figure 4: 



For r —* 00, after a certain number of oscillations, the pendulum falls down in the stable 
equilibrium point, see Fig.5, corresponding to a = (2p + l)ir, p integer, due to the force 
(1 + 2/ ) and to the damping term. This situation does not correspond to SBC, since 
a = (2p+ 1)tt for r — > 00. Rather, it gives a copy Ai obeying WBC. Indeed, from eq.([47|) 
we get 

~ ( M (cos(2p + l)vr) - 1\ 1 1 

~ const. + ^ n • (54) 

\ 2 J r r 



Figure 5: 
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• ii) The second case corresponds to a force / which, for a sufficiently large interval of 
time r, is negative enough, i.e., f < —1/2, see Fig. 6. In this case we have the following 
configuration, see Fig. 7. The pendulum starts at r — > —oo from an unstable position, 
a = 2irm, with velocity a' ~ 7e r . After some oscillations, it can come back to the 
unstable position, 2-7rn, at r — > oo, under the effect of the force /, and finally it can remain 
there, see Fig.8. This situation corresponds to SBC, implying that the field Ai decays 
faster than 1/r. 




Figure 6: 



a 




Figure 7: 
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a=2 K n 



Figure 8: 



• Summary 

This example shows that, starting from the configuration 

i 

Ai = ^eijkXjakfir) , diAi = , (55) 

we can obtain an equivalent field 

Ai = SidiS + &AiS , 9 4 1 4 = 0, (56) 

compatible with both WBC and SBC. 
In the case of SBC 

a(r) r ^oo — ► 2nn , S r -,oo — ► const. (57) 
rA] . 



For WBC 



a{r) r ^oo ->• (2p+ 1) 7T , 
<SV-»oo — ► cr ■ n , S depends on the orientation, S T -,oo — ► <5(#, </>) 

1™ ~ - • (58) 
r 
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2.2.1 The winding number of the Gribov copies 

Let us evaluate now the winding number, see Appendix C, of the gauge transformation S 

2 2 v ; 

We shall discuss in detail the case of SBC, namely 

a(r) r ^oo — * 2-nn , (60) 

so that 

SV-joo — > const. (61) 

Equation (|61|) implies that the space E? turns out to be compactified to the sphere S 3 , by 
identifying all points at infinity. We have to evaluate 

(62) 
(63) 
(64) 



From 



24^ Tr / dx ^ dx ^ (sfldiS) ^tys) (&d k s) . 
tfdiS = - (diSA s , 



we get 

v = -^Tr J d 3 xe ijk (ftS) (s^djS) (d k S^ 
It is useful to employ the notation 

S = N 4 + iN , S t = N 4 -iN -~a 

a ■ a 

N 4 = cos — , Ni = m sin — , 

N a N a = 1, a = 1,2,3,4. 



(65) 



In what follows we shall take in due account the surface terms originating from integration by 
parts. Thus 

1 



24vr 2 



Tr / d 6 xe ijk diN 4 + i diN ■a)) (N 4 -iN-a 



x ( [ d j N 4 + i(d j N ■a)) [d k N 4 - i (d k N ■ a 
1 



n -TV / d x£jA k 
24vr 2 ' J 



diN 4 [N 4 - iN ■ a + diN iN 4 + N ■ ~a 



Since 
we have 



x [djN A [d k N 4 - id k N ■ a + [djN ■ a) [id k N 4 + d k N ■ a 



e ijk djN 4 d k N 4 = , 



Tr 



(66) 
(67) 



diN 4 ^N 4 -iN -a) l^djN ■ a) ld k N ■ a 
+N 4 (diN • ^) ((d 3 N 4 ) (d k N • ^) - (d k N 4 ) (djN ■ 

'diN • ~a) (iN 4 (djN • ~a) (d k N • ~a ) - i {djN 4 ) (n ■ ~a ) ( d k N • ~a 



+ {diN--a) (N--a) (djN-^) [id k N 4 + d k N ■ a' 



(68) 
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From 

OiOj = 5ij + iEijk<Jk 

we find 



TroiOj = 25ij , (70) 
Tr (<Ti<Tj<r k ) = 2ie ijk , 
Tr ((ji(jjG k ap) = 2(5ij5 kp - 5 ik 5j P + 5 ip 5j k ) . 



Furthermore, observing that 



0~mO~pO~qO~r 



it follows 



1 



Tre i3k {8iN m ) N p (djN q ) (d k N r ) 
2e ijk (diN m ) N p {djN q ) (d k N r ) (5 mp 5 qr - 5 mq 5 pr + 5 mr 5pg ) 
2e ijk [(diN m ) N m (djN q ) (d k N q ) - (cW m ) N p (8jN m ) (d k N p )} 
+2e ijk (diN m ) N p (djN p ) (0 k N m ) 

, (71) 



£>xe ijk e m pq [2N A {diN m ) (djNp) (d k N q ) - 6 ($JV 4 ) N m (d 3 N p ) (d k N q )\ . (72) 



24tt 2 

Integrating by parts the second term, one obtains 



1 ' ,3, 



v = J drxeijhEmpqNi (diN m ) (djN p ) (d k N q ) 

~4^2 / ^ x£ ijkdi(e mpq N i N m (djNp)(d k Nq)) . (73) 
Concerning the surface term in ea. (|73|) . it turns out that 

£ijk£mpqNiN m (djNp) (d k N q ) 

= cos ^ sin %i m [dj (n p sin ^ ) [d k (n q sin ^ J ) e ijk e m pq 



a . 3 a 



cos — sin 



2 2 



flm (djflp) (d k flq) ^ijk^mpq 



1 ol . 3 a 

~2 cos — sin — n m (dj p - rijn p ) (d kq - n k n g ) Eij k £ mpq 
1 a o a 2 a o a 

^2 cos - sin -n m E ijk e m j k = — cos - sm —<m ■ (74) 



o f n f 2 Ot n CY 

(TxEij k di (EmpqN^Nm {djN p ) (d k N q )) = / <Fxdi I cos - sin" 1 -? 



ni . (75) 



Thus 



Moreover, recalling that a(r) r ^oo — > 2im, it follows that the surface term (|75|) vanishes. 
• Remark 

It is worth noticing that the surface term ([75 1) also vanishes in the case of WBC , due 
to the presence of cos^. Indeed, for WBC, a(r) r __ >00 —* (2p + l)7r, so that cos = 
cos(p7r + 5) = 0. Thus, for both SBC and WBC, we have 

1 r*„ 



v = ^2 / d 6 xE ijk e mpq N 4 (diNm) (djNp) (d k N q ) . (76) 
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Let us proceed then with the evaluation of expression (|7fi|). We have 



^2 J d 3 X£ij k£mpq cos ^ (<9* (n m sin (d,- 



n D sin ■ 



( ra sin 



a 



^2 / « xe ijk e mpq cos - ^sin - (d;n m ) + — cos - (d,a 



n 



7) 



x (sin - (djUp) + cos - (5,-a) ) (sin - (d k n q ) + cos - (<9 fc a 



n 



n,, 



o 



1 

3^ 



(f X£ij k £ m pq COS 



2 L 



sin 



(diUm) (djUp) + 



a 



a 



sin ■ 



cos - (<9jn m ) (9j-a) 



a; ck 



^2 y d 3 xEij k £mp q cos ^ (sin 3 ^ (<9jn TO ) (djiip) (d k n q ) 



+ sin 2 - cos ^ y (<9;n m ) (<9jn p ) (5 fc a) 



sin 2 ^ cos ^ (din m ) (dja) (d k n q ) 
2 a 



sm z - cos - (dja:) (<9,n p ) (d&n g ) 

^2 y d 3 x£i jk £ mpq cos ^ (sin 3 ^ (<9jn TO ) (<9,n p ) ((9 fc n g ) 

3 . 2 a a WQ w _ , 

+ 2 sin - cos -n q (din m ) (djUp) (d k a) 



Furthermore 



y d3x ~3 cos f sin3 ~^ £ ijk £ mpq ($im ~ riin m ) (5 jp - UjUp) (5 kq - n k n q ) 

l/' I q.2 a 2 a /l I c w r \ 

+ 27r2 / "2 2" ~2 £i i k£m P1 n l nk (° im ~ n i n ™> (°3P ~ n j n p) 

2 / d X—^ COS — sin ~£ijk^mpq {.(&im$jp ^im^j^p ^jpf^i^m) (&kq ^k^q)) 



3vH J r 6 2 2 

1 /" n . O nfl 

H ^ / a x sin — cos - ^ 

2vr 2 ./ 2 2 r 

1 

3^2 



1 ^ r/ 3 ■ 2 — 2 — '— (XX _ A _A ^ 

ffl ISin —COS 77*-^ ~2~£'ijk£mpq'H i q 'n k [Oi m 0jp <Ji m Tljnp Oj p Tlin m ) 



d X „ COS sin „£ijk^mpq (&im i&jp&kq ^jpf^k^q ^kqf^j^p) ^jp^kq^i^r, 



+ 2~^2 / d xsm 2~" COS 2"° ~2 £ ijk £ mpqn q n k d im dj p 

1 /" 3 1 a . 3 a (xxx "ix X \ 

2^_2 / "^^3 "2" ~2~£ijk£mpq {Oi m 0jp0 kq oOi m OjpTl k nq) 

1 /" ,0 . 2 a 2 a /l 

+ 27r2 / m "2 COS 2 a ^ £ijk£ijqnqnh 

1 /" , 1 a , o a 

^2 / « cos 2 sm "2" y £i i k£i i k ~ ?>£ijk£ijqn k nq) 



1 r ,0 . 2 a 1 a 1 1 
H — 5- d xsm — cos —a 
7r 7 2 2 r z 

1 /" j3 ■ 2 « 2 « / 1 
a xsm —cos —a -77 
7r^ / 2 2 r z 
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In order to evaluate this last integral we shall make use of polar coordinates 

v = — / dr sir — cos —a = / dra e 2 + e 2 e 2 - e 2 

Wo 2 2 47r7 V J \ J 

= -— dra' (e ta + e~ ia + 2) (e ta + e~ ia - 2) 

1 Z" 00 

= — / dra' (1- cos 2a) . (79) 
2?r Jo 

Finally, for the winding number v associated to the gauge transformation S, we find 

1 1 

v = — (a(oo) - a(0)) (sin2a(oo) - sin2a(0)) . (80) 

27T 47T 



Summary 

Let us consider first the case of SBC. From 

a(r) r ^oo — > 27r?i , n integer 

a(r) r ^o — > 27rm , m integer , (81) 

follows 

vsbc = integer (82) 



For W^BC 



a(r) r ^oo — ► (2p + 1) 7T , p integer 

a(r) r ^Q — > 2-7rm , m integer , (83) 

^BC = (integer + ^ . (84) 



• Remark 

It should be noted that, in the case of WBC, the expression Q6'#j) does not have the meaning 
of a winding number. In fact, due to the behavior of the gauge transformation S at infinity, 
eo.([5flJ). the space R 3 cannot be compactified to the sphere S 3 . As a consequence, expression 
(I6'ifl may now take non integer values. 
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2.2.2 The Gribov copies of the vacuum 



Let us discuss now the existence of the Gribov copies of the vacuum, Ai = 0, which obviously 
satisfies diAi = 0. We look at the equivalent field 

7 „tan 1 ■ 9n fl l\^dn 1 , 

Ai = b'Oib =-sma— M-eosa ra— 1 — a (r)nm , 

2 dxi \2 2 J dxi 2 w % ' 

5 = e H } ')"-^ = cos|+^-"^sm|, (85) 
which corresponds to a pure gauge 

A i = tfd i S, F ij (A) = 0. 
The gauge condition diAi = reads now 

2 2 

a (r) H — a (r) ^ sin a = , (86) 

or, r = logr, 

corresponding to the damped pendulum of Fig. 9. 




Figure 9: 



The regularity condition of the equivalent field Ai at the origin, r = 0, r = — do, gives 

a(r) — ► 27rm + 7r /or r — ► , m integer (88) 
a(r) — > 2irm + 7e T /or r — > — oo , 

which holds for both WBC and SBC. Since the force / is now absent, the pendulum starts at 
r — > — oo from the unstable position, a = 2nm, with velocity a 1 = r ye T , see Fig. 10. 
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Figure 10: 



After a certain number of oscillations it comes to the stable position a T - 
the constant force and to the damping, see Fig. 11. 



(2p + 1) 7r, due to 



!-}+ -r. 



Figure 11: 



Thus, in the absence of the force /, there are no Gribov copies of the vacuum if SBC are 
imposed. However, if WBC are adopted, even the vacuum has Gribov copies. As we have 
already seen, this corresponds to an equivalent field Ai which, for a T -,oo = (2p + 1) vr, behaves 
as 

Ar^oo ~ (cos(2p + 1)tt - 1) - ~ - . (89) 



19 



2.3 The Henyey example 

In this section we shall discuss a class of Gribov's copies proposed by Henyey 0. The relevance 
of Henyey's work is due to the fact that Gribov's copies with vanishing winding number and 
which fall off faster than 1/r, for r — ► oo, are explicitly obtained. 

The starting configuration is 

~A = i~aa 3 , d, L Ai = . (90) 

Notice that the gauge configuration, ea.(|90|). lies in the abelian diagonal subgroup of SU(2). We 
shall consider the set of gauge transformations parametrized by 



S = e iaf ^ = cosa + iJ -~asina , (91) 



where / is the unit vector 



/ • / = 1 , (92) 
/ = ~e x cos/3 + "e^sin/3 , 

and (~e x , ~~Zy> ~^z) are the unit orthogonal cartesian vectors. Introducing the vector ~~g 

g = — e x sin (5 + ~~e y cos j3 , (93) 



it turns out that the set ( / , g , e z j yields a right-handed orthonormal triad which rotates 
about the z-axis. Also, 

dil = -gdi(3 , 8^ = -7di(3 . (94) 
Let us evaluate the gauge transformed field A{ 

M = S^diS + S j AiS (95) 
= i ^cos a — i f -~a sin aj aio 3 ^cos a + i f -~a sin aj 

+ ^cos a — i f -~o sin aj (— sin a + i f -~a cos aj dia 

+i (cos a — i f -~a sin a ) ~~g • ~o sin adi/3 . 



Recalling that ( / • ~a ) = 1, it follows that 



Ai = ia,i (cos 2 aa 3 + i sin a cos af k [a 3 , a k ] + sin 2 afkfmCTkVsVm) (96) 
+id. t a f -~a + di[5 (i sin a cos a~g ■ ~~a + sin 2 a ( f ■ (~g ■ 
= iai (cos 2 aa 3 - 2 sin a cos ae 3km a m f k + sin 2 af k f m a k (5 3m + ie 3mp a p )) 
+idia f -~a +difi(i sin a cos a~g ■ ~a + sin 2 afig k (5 ik + ie ikm a m )) . 
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From 



£3kmV m fk = ( e z x f j ■ a = (e z x (e x cos /3 + e y sin /?) ) • cr 
= (e y cos j3 — e x sin /?) - a = g - a , 
figk£ikmVm = (f x ~g) ■ ~o = ((~e x cos (3 + ~e y sin/3) x (-~e x sin(3 + ~e y cos/?)) • ~a 



= e 2 • a = a 3 , 

f k f m£'Arnp@ k@ p — fkfm^3mp ($kp ~i~ i^kpq&q) — fkfm^3mk ^fkfm^Smp^pkq^q 
= —ifkfm(53k5mq — 83 q 5 m k) CTq = i(J3 y 

h = 0, 7-^ = 0, (97) 

it follows 

Ai = icii (cos 2 ao3 — 2 sin a cos a^f • 7? — sin 2 aas) + ic^a / • ~~o 

+dif3 (i sin a cos a; if • ~~o + ? sin 2 ao^) . (98) 

Finally, from (cos 2 a — sin 2 a) = cos 2a, we get 

= i(T3 (aj cos 2a + 9j/?sin 2 a) + «5ja / ■ ~a 

+i~~g ■ ~~o ai sin 2a + -9j/3sin2a^ . (99) 

Let us turn now to the condition diAi = 0. Computing the divergence of Ai, yields 

diAi = ia 3 (— 2aidia sin 2a + d 2 p sin 2 a + 2 sin a cos adiPdia) 
+id 2 a f • ~a + idiadiP~~g ■ ~a 
—idiP f ■ 7j ( —a-i sin 2a + - dip sin 2a 



-/ // • t j —2aidiCt cos 2a + -^p 2 P sin 2a + di(3dia cos 2a) 



Following Henyey [2], we shall impose 



I 



(100) 



d 2 (3 = didi/3 = , aidia = , d^a = . (101) 
Thus, the condition diAi = reduces to 

d 2 a + aidiP sin 2a - ^diPdiP sin 2a = . (102) 



These conditions are more restrictive than necessary. 
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• Summary 

A = i a (73 , diAi = diOi = . (103) 
S = e ia ' ' u = cos a + i f -~a sin a , 

7 = "e !E cos/? + "ej / sin/? , f.f=l. (104) 

-For the gauge transformed field we get 

Ai = S^diS + S^AiS 

= io~3 (a{ cos 2a + di(3 sin 2 a) + idiot f ■ ~a 

+i~~g ■ ~o {^— ai sin 2a + —di(3sm2aj . (105) 

The condition for the existence of Gribov's copies, diAi = gives 

8 2 f3 = didip = , tudio = , dipdia = , (106) 

and 

d 2 a + aidip sin 2a - -di(3di/3 sin 2a = . (107) 

2.3.1 Henyey's solution 

In order to solve the eqs. ljl06J ). 1)107(1 it is useful to adopt polar coordinates |2}, see Appendix D. 
We set 

~a = a(r, 0)~i% , j3 = /%) , a = a(r, (9) . (108) 

The equations 



V-~a = —^-^-a(r,9) = 0, (109) 
r sin dip 

_> -> i a 

a • Vet = a(r, 0) -— — 0) = , 

r sin cV 

Va • V/5 = , 



are fulfilled. From 



V 2 /3 = -r^^/3 = 0, (HO) 



9 ■ 2 



0d^ : 



it follows that we may take 

/3 = <p. (Ill) 
It remains to solve the condition ()1U7(I . which now reads 

1 3 / 2 3a\ l d ( . „da\ . _ ( a \ 



ri Or \ r fr) + 50 {^W J + ^ [T^O " rt^ 1 = ° " ' 

Setting 

a(r, 0) = rfe(r) sin , (113) 
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equation (j!12j) becomes 

b + (r 2 b" + Arb') sin 2 9 + (a - ^T^g ) sin ( 2r6 sin = ■ ( 114 ) 

The strategy adopted by Henyey [2] in order to solve this equation is that of expressing o in terms 
of b, and then searching for a suitable b which gives the desired behavior for a. Accordingly, we 
write 



a 



2r sin 9 sin (2rb sin I 
We look then for a function b such that: 



(b + r 2 sin 2 9 (b" + 4-U . (115) 



a is regular at the origin, r — ► , (116) 
a decays faster than 1/r for r — > oo . 

First of all we notice that expression (j!15|) is not singular at = 0. Although each term of 
eq, (|115[) is singular, their sum has no singularity at 9 = 0, i.e. 

Also, in order to avoid possible singularities in the term (1/sin (2r6sin#)) we require 

2rb(r) < vr . (118) 

Since sin# < 1, the condition (|118|) implies that the argument of sin(2r&sin#) cannot be equal 
to 7T, i.e. 2rbs'm9 < ir. 

Let us look now at r = 0. We have 

1 1 (, . o . A„ .b' 



2rsin# 2r6sin# 
Requiring that 



^6 + r 2 sin 2 (&" + 4^U . (119) 



b(r) r ^o ~ b o + hr 2 + 0{r A ) , with b ^ , (120) 

it follows 

Or-+0 « ^ p 7T " + °W « 0(r) , (121) 
2r sin # zro smt^ 

so that a is regular at the origin. It remains now to discuss the limit r — > oo. In this case we 
search for a &(r) such that r 2 a — > for r — ► oo. Let us set 

( 122 ) 



and let us determine n. From 

„2„ ^ „,2 



1 1 1 \ r 4 sin 2 6l 



r a r ^oo ~ r /n . „. .„ . „. (n(n + 1) — 4n)) — 

\2rsm9 r n sm(^°)J sin(^)^ 7 " r' 

2 / 1 1 1 \ , „ s r 4 sin 2 ' 



Thus, if n = 3, 



r (^2rsin0 r«(M^)+0(^)J n( " 3) r™+ 2 sin (2£#) 
„ / 1 \ . . r sin 

°(^)-^- 3 )— • ( 123 ) 



&(0r-oo « ^> ( 124 ) 



23 



Summary 

In summary, any function b(r) such that 



b(r) r ^ ~ b + b 2 r 2 + 0(r 3 ) , with b ^ , 
2rb(r) < tt , 

KO™ ~ 3, (125) 



r 3 

will give a gauge field A 

A = ill '03 , diAi = diai = . (126) 
~a = a(r, 0)~e v 

air, 9) = l — -4 777 (b + r 2 sin 2 9 ( b" + 4- 

2r sinf sm(2rosm#) \ \ r 

which is regular at the origin, r = 0, and decays faster than 1/r 2 at infinity, r — > oo 

a r _> ~ O(r) , 

^ 2 «(0™ « ofzV (127) 



An example of such a function b(r) is given by 

b(r) = 



(r 2 + r^) 



2\3/2 



* < T i (128) 
where the condition k < (3 3 / 2 7rr 2 ) /4 stems from the requirement 2rb(r) < tt. 

2.3.2 The winding number of Henyey's solution 

It remains now to discuss the winding number of Henyey's solution. Let us begin with the 
evaluation of the winding number corresponding to the gauge transformation ()91|) 

S = e ia ~f^ = cos a + if ■ ~a sin a . (129) 
From expression (|73[). we have 

/ d 6 xe ijk e mpq N A (diN m ) (djN p ) (d k N q ) 

d*xe ijk d t (e mpq N 4 N m (djNp) (d k N q )) , (130) 



1_ / ,>> 

3tt 2 . 

1 

"4^2 



with 



Recalling that 



N4 = cos a , Ni = fi sin a . (131) 
dif = ~gd t p , 7-^ = 0, (132) 
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it follows 

v = ^ j d3x£ ijk£m Pq cos a [(sin ag m dif3 + / m cos a^a) (sin ag p dj[3 + / p cos adja) 
x (sin ag q d k (3 + / 9 cos a<9 fc a)] 
~4^2 / d 3 x£ijkdi[E m pq sin a cos a/ m (sin ag p dj/3 + / p cos adja) 

x (sin ag q d k j3 + / g cos ac^a)] (133) 
= / <^ x£ ijk £ m P q cos a [(sin a cos ag m f p di{3dja + sin a cos af m g p diadj(3) 

x (sinagg<9 fc /3 + / g cos a<9 fc a)] 
~4^2 y d 3 xeijkdi(e mpq sin a cos af m (sin ag p dj/3) (sin ag q d k f3)) 
= . (134) 

We proceed now with the computation of the winding number of the starting gauge configuration 
(|9"T|) and of the gauge copy (|95|). To obtain the winding number of the expression (|§T|) . we start 
from the Pontryagin index, eas. ([376|) — (|378f) of Appendix C, 



„-Tr / dS ^upcf I Aj,<9p^4 CT + —A u ApA a ) . (135) 

87r -/S3, V 3 y 

From 

~A = i~aa 3 , diA { = , (136) 

and 

A 4 = , (137) 

expression (|135|) becomes 

v = ~^2 Tr J s3 dS A e Uj k yAidjA k + ^AiAjA^ . (138) 

Since the starting field is in the abelian subgroup of SU(2), it follows that 

TrAiAjA k = -iaiaja k Tr (cr 3 a 3 a 3 ) = . (139) 

Also _^ 

£Mjk T r (AidjA k ) = -2e Mjk aidja k = -2~a ■ ( V x ~a ) . (140) 

Furthermore, since ~~a = a(r,9)~e ip , eq. ([126|) . has component only along the direction e v and 
V x ~~a has components along the directions e r , ~eg, it follows that 

a^Vx a)=0. (141) 

We see therefore that the Pontryagin density £aj k (AidjA k + ^AiAjA k ) associated to the start- 
ing gauge configuration vanishes. Finally, for the winding number of the copy (|95|) one has 

^copy — ^original ~i~ ^transformation — . (142) 

Therefore, we have shown that it is possible to construct Gribov copies with vanishing winding 
number and which fall off faster than 1/r for r — > oo. This concludes the discussion of Henyey's 
example. 



25 



3 Part II: The Gribov Horizons 



3.1 Generalities 

In order to introduce the notion of Gribov horizon let us look at the eigenvalues of the Faddeev- 
Popov operator, i.e. 

-d^(d^+[A^]) = e(A)^ . (143) 

As already underlined in Part I, equation l)143|) can be seen as a kind of Schrodinger equation, 
with playing the role of the potential. For small values of A^, ea. (|143|) is solvable for positive 
e only. More precisely, denoting by e%(A), 62(A), 63(A), the eigenvalues corresponding to a 
given field configuration A^, one has that, for small A„, all €i(A) are positive, i.e. €i(A) > 0. 
However, for a sufficiently large value of the field A^, one of the eigenvalues, say ei(A), turns 
out to vanish, becoming negative as the field increases further. As in the case of the Schrodinger 
equation, this means that the field A^ is large enough to ensure the existence of negative energy 
solutions, i.e. bound states. For a greater magnitude of the field A^, a second eigenvalue, say 
62(A), will vanish, becoming negative as the field increases again. Following Gribov [Q, we 
may thus divide the functional space of the fields into regions Co, C\, C2, C n over which the 
Faddeev-Popov operator, — <9 M (d^ ■ + [A^, •]), has 0, 1, 2, n negative eigenvalues, see Fig.12. 




Figure 12: The Gribov horizons 



These regions are separated by lines h, I2, h, l n on which the Faddeev-Popov operator has 
zero energy solutions. The meaning of Fig.12 is as follows. In the region Cq all eigenvalues of 
the Faddeev-Popov operator are positive, i.e. —d^ (d^ • + [A^, •]) > 0. At the boundary l\ of 
the region Cq the first vanishing eigenvalue appears, namely on l\ the Faddeev-Popov operator 
possesses a normalizable zero mode x 

d fl (d llX +[A fl , X })=0. (144) 

In the region C\ the Faddeev-Popov operator has one bound state, i.e. one negative energy 
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solution. At the boundary I2, a zero eigenvalue reappears. In the region C2 the Faddeev- Popov 
operator has two bound states, i.e. two negative energy solutions. On I3 a zero eigenvalue shows 
up again, and so on. The boundaries l\, 1%, I3, l n , on which the Faddeev- Popov operator 
has zero eigenvalues are called Gribov horizons. In particular, the boundary l\ where the first 
vanishing eigenvalue appears is called the first horizon. 

• Remark 

It is useful to emphasize that in the region C$ the Faddeev- Popov operator has only positive 
eigenvalues. Therefore, this region can be defined as the set of all transverse fields for which 
the Faddeev-Popov operator is positive definite, namely 

C = {Ap, 8A = 0, -dp (dp • + [Ap, •]) > 0} . 

• Remark 

In order to obtain a better understanding of the notion of the Gribov horizons, let us remark 
that there is a close relationship between the horizons and the existence of Gribov copies. 
In part I, we have discussed the existence of equivalent fields by considering finite gauge 
transformations 

Ap = S^dpS + rfApS . (145) 
The requirement that the field Ap obeys the same transversality condition as Ap 

dA = dA = 0, (146) 

yields the equation 

d^dpS + s^d^s + b^a^s+^a^s = o . (147) 

For S close to unit, S = 1 + a, aCl, expression (j!47j) reduces to 

dp(d fl a + [A f ,,a]) = 0. (148) 

We see therefore that the condition for the existence of an equivalent field Ap close to Ap, 
i.e. 

Ap = Ap + (d^a + [Ap, a]) , (149) 
relies on the existence of a zero mode for the Faddeev-Popov operator. 

• Remark 

The transversality condition of the Landau gauge, dpAp = 0, implies that the space-time 
derivative dp and the covariant derivative Dp(A) = (dp ■ + [Ap, •]) obey the following com- 
mutation relation 

dpDp(A) = Dp(A)dp . (150) 

As a consequence, the Faddeev-Popov operator M = —dpDp(A) turns out to be Hermitian, 
M = Jvft . Its eigenvalues are thus real. In fact, integrating by parts, one has 



d 4 x ((dpDp^ <p a = I d*x (dpdp^ + f abc A b ^dp^) <f a 

d 4 x (r ] dpdp<p a - f abc ^Aldp<p a ] 
d^x ^dpdpip a + f abc r ] 'A%<p c 
j atx^ (dpDpip a ) . (151) 
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3.2 Example of a zero mode of the Faddeev-Popov operator 

It is useful to provide here an example of a normalizable zero mode of the Faddeev-Popov 
operator. We shall work in three dimensions, the gauge group being SU(2). The aim is to 
obtain a normalizable zero mode x-> solution of 

V 2 X + [A i ,d lX }=0. (152) 

We shall follow Henyey's strategy [2], see also and write 

~A = i~aa 3 , d { Ai = . (153) 

Adopting polar coordinates, see App.D, we shall set 

~a = a(r, 9)~e ^ . (154) 

We look now at a zero mode x °f the form 

X = a(r, 9) (<7i cos p + 02 sin (p) . (155) 



From 
it follows 

Therefore 



Setting 
we get 



[<Ti,(Tj] = 2iEij k (Tk , (156) 

[0-1,02] = 2ia 3 , [0-3,0-1] = 2io 2 , [o 3 ,o 2 ] = -2ioi . (157) 

[Ai,dix] = i [a 3 ,ai]aidi (a cos f) + i [a 3 ,a 2 ]aidi (a sin if) 
= —2o2aidi (a cos if) + 2o\aid{ (a sin ip) 

a d (a cos ip) a d (a sin ip) 

= -2o 2 — — ^ h2oi — — . 

rsmtt op rsmO op 

(158) 

a(r, 9) = rb(r) sin# , (159) 
[^4j,0«x] = 2absmpo2 + 2ab cos pa\ . (160) 
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Also, 



^2 19 / 2 dx\ 1 9 ( . n dx\ 1 d 2 x 



r 2 dr \ dr J r 2 sin 9 89 \ 89 J r 2 sin 2 9 dip 2 

■ n, ■ ^ 1 8 ( 2d(rb(r)) 

sin 9 [ax cos 99 + 02 sin ip) — I r 



r 2 <9r y 9r 

6(r) 8 ( 8 sin 6* 
+ (cti cos 99 + 02 sin 99) — ^ ™ sin" 



r sin 9 89 \ 89 
b(r) 8 2 (o\ cos ip + o"2 sin 99) 
r sin dtp 2 

1_ 8_ 

b 8 



sin (<7i cos 99 + 02 sin 99) — — (r 3 6' + r 2 6) 



+ (oi cos 93 + 0-2 sin 99) — : — - — (sin 9 cos I 
r sin ft at> 

[a 1 cos 99 + (J2 sm 99) 



r sm ( 

sin 6* (0-1 cos ip + a 2 sin 99) — (4r 2 b' + r 3 6" + 2rb) 

+ (a 1 cos 9? + 02 sin 99) (cos 2 9 — sin 2 9) 

rsmv 

b{r) 



(a 1 cos 99 + (72 sin 99) 



r sm ( 

= (46' + r&") sin 9 [a\ cos 99 + 02 sin 99) 

6 / (cos 2 6 - sin 2 9) \ 

+ - 2 sin 9 -\ (<ti cos 99 + o"2 sin 99) 

r \ suit) / 

, • V 

: — t (cti cos 99 + cr 2 sm 99) 

r sm 

= (46' + rfr ) sin 9 (a± cos 99 + 172 sin 99) 

b b(r) 

H ; — r ((7i cos 99 + cr 2 sin 99) — - (<7i cos 99 + cr 2 sin 99) 

r sm r sm # 

= (46' + r6") sin 6» ((7i cos 99 + <7 2 sin 99) (161) 

Finally, 

V 2 x = (46' + rb") sin (crj cos 99 + <7 2 sin 99) . (162) 
The equation H152|) is thus equivalent to 

(46' + rb") sin 9 (01 cos 99 + 172 sin 99) + 2ab sin 99172 + lab cos 99171 = , (163) 

which yields 

(4b' + rb") sin 9 + 2ab = . (164) 
Following Henyey |2], for a we obtain 

a = sin e(6" + ^) . (165) 

We now look for a function 6 which yields a field configuration a(r, 0) which is regular at the 
origin, r = 0, and which decays faster than 1/r for r — > 00. As done in Part I, for b we choose 

k 

b(r) = -7- , k const, and rn 7^ . (166) 

(r 2 + r 2 ) 3/2 
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Therefore, for a(r, 9) we find 6 



/ /\ \ r , g o\3/2 „ / d r 12 

a(r,6») = --(r 2 + rl) ' sin (9 -3 



2 7 " V 8r {r 2 + rlf' 2 (r 2 + r 2 ) 

i:( r 2 + r 2) 3 / 2 siri ^ 1 5. /J 



lr 2 



+ r 2)7/2 ( r 2 + r 2)V2^ 



2 



- rr ° 2 sinfl. (167) 

2 _|_ ^.2^ 



Observe that 



a(r, 0) -> for r , (168) 
a(r, 9) ~ -r for r — > oo . 

It remains to check that x ls a normalizable zero mode with respect to the scalar product 

(X | Y) = Tr j d 3 xX j Y . (169) 

From 

X = rb sin 9 {a\ cos (p + a% sin 93) , (170) 

we have 

t^t = rb sin (o"i cos tp + 0% sin </>) , (171) 

and 



Tr 



(x^xj = 2r 2 b 2 sin 2 9 . (172) 
Thus 

//■7T (>oo j.4 
d 3 xr 2 b 2 sin 2 9 = 4nk 2 sin 3 9 d9 dr 3 < 00 , (173) 
Jo Jo (r 2 + r 2 ) 

showing that we have found a normalizable zero mode of the Faddeev-Popov operator. 
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• Summary 

For the zero mode we have 

V 2 x + [A i ,d iX }=0, (174) 



with 



and 



A = i~aa 3 , diAi = , 

~a = a(r,0)~e,p , (175) 

X = rb sin 9 (<7i cos ip + 02 sin ip) , (176) 
k 

b(r) = , k const, and r / . (177) 

(r 2 + rl) 

(x\x) = J d 3 xTr ( X f x) < °o • (178) 



For a(r, 6) one finds 



26 V r J 2 ( r 2 + r 2) 2 

a(r,0) -► for r -»• , (180) 
a(r, 0) ~ — for r — ► oo . 

Also, a is transverse 

1 da(r,9) , N 

9iOi= . a \ = • 181 
r sin a<£ 
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3.3 An important statement 



Let us prove the following statement PP: 



• Statement 

For any field located within the region Ci, close to the boundary l\, there is an equivalent 
field within the region Co, close to the same boundary li, see Fig. 13 











^^^^^ 













Figure 13: The Equivalent Fields 



Let be the field located in Co, close to the first horizon l\. We write 

A/_ t = C M + , a M small with respect to C M , (182) 

where the transverse field C^, dC = 0, lies on the Gribov horizon l±, i.e. it exists a normalizable 
zero mode (fo such that 

^(^o + [C M ,^o]) = 0. (183) 

• Remark 

Notice also that, by definition, the field A^ is located in the region Cq, so that 

-d lt (a it - + [A ttr ])>0. (184) 

The eigenvalues problem for the Faddeev-Popov operator corresponding to the field A^ can be 
easily handled by means of the perturbation theory of quantum mechanics. Let us look indeed 
at the eigenvalues equation 

-Bp (d^ + [C„ + o M , ij>]) = e(a)i/> , (185) 

i.e. 

- {d 2 ^ + [C M , d^} + [a,*, 1>]) = e(a)V . (186) 

One might think of as a small perturbation. Accordingly, the shift e(o) of the eigenvalue of 
the Faddeev-Popov operator from zero can be obtained in perturbation theory by evaluating the 
expectation value of the term [a M , •], namely 

, x (<Po I fy. [a^ipo]) Tr J d^xiifodf, [a^ipo]) 

61(21 — ; j : — — =— j-j 7 s • I J-0( ] 

(ip I (po) Tr J d 4 x (vwo) 
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Since 9 M (d^ifo + [C M , ipo\) = 0, it is easily verified that the field 

A IM = A li + D lt (C)<p , (188) 

has the same divergence of A^. In fact 

= + d^D^C)ip = . (189) 

The fields A^ and A^ can thus define equivalent fields. For that, we have to find a gauge 
transformation S such that 

A„ = St^fif + rfApS . (190) 
We shall look at S close to unit, in the form 

S = 1 + a + ^a 2 + .... (191) 



Also 



Indeed 



5t = l- a + I a 2 + .... (192) 



= 1 - a + ^a 2 + a - a 2 + ^a 2 + 0(a 3 ) 

= l + 0(a 3 ). (193) 



For A M we get 



A u = [l-a + la 2 \ A u [l + a + ^a 2 ) + [l-a + ^a 2 )d u [l + a + ^a' 



From 
it follows 



2 

= (l-a + ^a 2 ^j (a^ + A^a + ^A^oA 

+ ^1 - a + ^a 2 ^ ^<9 M a + ^ (d^a) a + (d^a) 

1 1 

= Afj, + v4 M a + -A^a 2 - aA M - ai^a + -a 2 A^ + d^a 

11 o 

= A M + -D M (A)a + i (<9 M a) a - (d^a) 

+^A fI a 2 - aA^a + -a 2 ^ + 0(a 3 ) . (194) 



[a, <9 M a + L4 M , a]] = a (d^a) - (d^a) a + aA^a - a 2 A^ - A^a 2 + aA^a , (195) 
\ = A IX + D fl (A)a - X - [a, d^a + [A^ a}} + 0(a 3 ) . (196) 
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• Remark 

We must retain the second order term in the equation (|191|) otherwise, from 

S = 1 + a , (197) 

we would obtain 

A lx = A„ + D„(A)a. (198) 
Condition (|189|) . dA = dA, would thus imply 

d^D^Aja = , (199) 

which has no solution for a since A^ is not located on a horizon, see eg. (|184|) . 

From equation (|196|) one obviously has 

D^C)^ = D^{A)a - 1 [a, 8^a + [Ap, a]] + O(o 3 ) . (200) 

Tacking the divergence of both sides of ea. (j200|) . we obtain the condition to be fulfilled in order 
that A/j, and A^ have the same divergence, dA = dA, i.e. 

8MA)a - -0„ [a, d^a + [A^, a]] = . (201) 
This condition can be analyzed iteratively by setting 

a = (£>o + 5 , a small with respect to ipo , 

A^ = C M + , a M small with respect to A^ , (202) 

It is useful to introduce an expansion parameter A which we shall set to one at the end§ 

a = Xifo + A 2 5 , 
Ay = C^ + Xa^. (203) 



Condition (|2U1|) gives 

d 2 (Av?o + A 2 5) + d )Jl [C lx + Aa M , X^ Q + A 2 5] = [X<p + A 2 5 , (Av?o + A 2 5)] 

[XifQ + A 2 5 , [C M + Aa /t , A^ + A 2 5] ] 

(204) 

so that, up to terms of the order A 4 

A 2 d /t (d^a + [ C M , 5]) + A 2 d M [ a^, <p ] + A 3 <9^ [ o M , 5] = 4^ [</? Q , + [C M , 990]] 
+1^ bo, ^5 + [C M , 5]] + 4^ [5, a^^o + [C M , v?o]] + [^o, [o/x, </>o]] + 0(A 4 ) . 

(205) 



In particular, to the first order A 2 , we find 



d^D^C)a + [ o M , ^0] = ^ [^0, £> M (C>o] , (206) 



§The introduction of A turns out to be useful in order to analyze the condition (I20H . Observe indeed that, if 
A < 1, then A 2 < A. This means that a M is smaller than C M , and that a is smaller than ifio. 
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from which it follows that 



Tr J d A x (^D^C)a + ^ [ o M , <po]) = X -Tt f d 4 x [ip , D^C)y }) . (207) 

Moreover, due to property d^D^C) = 0^(0)8^, we have 

Tr f d A x m d^D^(C)a = Tr f d A x (d^D^Qtpo) 5 = 0. (208) 
As a consequence, condition (|207|) reads 

Tr J d 4 x (ipodp [ a„, <p ]) = X -Tr J d 4 x (ip 8^ fa, D^C)^)) . (209) 

• Summary 

The fields and An 

A ti = A lx + D„(C)<po, (210) 

are equivalent fields 

8^ = 8^ + d^D^C)^ = 8^ . (211) 
For the gauge transformation S 

A„ = S^d^S + SU^S , (212) 

we have 

S = l + a + ^a 2 + 0{a 3 ) , (213) 

so that 

Afj, = An + D )1 (A)a> - - [a, 8„a + [A M , a}} + 0(a 3 ) . (214) 
The condition dA = dA gives 

d^D^Aya = ^ [a, 8^a + a]] + 0(a 3 ) . (215) 



Setting 



a = (po + 5 , a small with respect to y?o ; 

-*V = C M + a M , a M small with respect to A^ , (216) 



we obtain, to the first order, 



d^D^a + <9 M [ a M , <po] = ^ fa, ^(C)^o] ■ (217) 



It remains now to check on which side of the horizon l\ the equivalent field A^ lies. Let us 
A„ = A„ + ZVCVo = C M + a M + £> M (C>o = C„ + < . (218) 
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rewrite as 



As done before, we evaluate the shift e(a') of the eigenvalue of the Faddeev-Popov operator from 
zero. Treating the field a'= + D^(C)(fo as a perturbation, one obtains 

/ /\ (go I d n Tr f rf 4 x (y a M [q^ yg] + ypg^ [D M (g)y> , y ]) 

CIO I — : j : — — 7. — —. ~, r . (Ziyi 

{cpo I ipo) Tr J d^x {ipocpo) 

Furthermore, from ea. (j2U9j) . it follows 

. ,. Tr f d 4 x (tpndu |a u , ipn]) , „ , 

Thus, if close to l\, is located in Co, e(a) > 0, there is an equivalent field, = A^+D^(C)ipo, 
close to Zi, which is located in C\, e(a') = — e(a) < 0. Also, it is worth mentioning that this 
derivation can be generalized to fields close to any horizon l n . This concludes the proof of the 
statement. 



3.4 Restriction of the domain of integration to the first horizon 

The previous statement suggests that the domain of integration in the path integral should be 
restricted to the first horizon, i.e. to the region Co where the Faddeev-Popov operator is positive 
definite. 

• Remark 

More precisely, two additional requirements should be fulfilled to justify the hypothesis that 
the domain of integration should be restricted to the region Cq. The first one is that one 
should be able to prove that not only for small neighborhoods close to the horizons l n , but 
also for any field in the region C n there is an equivalent field in the region C n -\. This 
would ensure that it is always possible to find a chain of gauge transformations which 
brings a field in the region C n to the corresponding equivalent field in the region Co . The 
second requirement is that one has to be sure that the region Co is free from Gribov copies. 
We shall come back on these issues later on in the Conclusion. For the time being, we 
shall assume that the significant range of integration in the path-integral is determined by 
the region Cq. 

Accordingly, for the partition function Z we write 



Z = M J DA fl DcDc5(dA)e-( SYM+ ^ d4x " ad ^ bcb )v(C ) 

= N f DA ll 5(dA)e~ SYM det [-d^5 ab - f abc A c ^ V(C ) , (221) 

where the presence of the factor V(Cq) means that the integration is performed only over the 
region Co- In order to characterize the quantity V(Cq), we look at the connected two-point ghost 
function (c a (x)c b (y)) . We have 



^{x)c b (y) 



N f DA^DcDc5(dA)c a (x)c b {y)e-( SYM+ ^ d4x ^W cb )v(C ) 



N I DA^5(dA)e~ SvM det (-8^) (d^D^) 



-i 



ab 



xy 



V(C ) . 



(222) 
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The presence of the factor V(Cq) in eq. ([222j) implies that (c a (x)c fc (y)) c can become large only 
when approaching the horizon l\. On the other hand, the behavior of the two-point ghost 
function obtained from perturbation theory, i.e. with V(Cb) = 1, is given by 

(lf{x)c b (y)) c = 6 ab f ^Q(ky^-v) 
1 



G(k) = nr, (223) 



where A is the ultraviolet cutoff and N is the Casimir of the adjoint representation of the gauge 
group SU(N) 

jacdjbcd _ N§ ab _ ( 224 ) 

The expression of G{k) in ea. (|223|l displays two singularities, at k 2 = and at k 2 = A 2 exp ^— p-^T 

However, the presence of the factor V(Co) makes it impossible for a singularity of G(k) to ex- 
ist at nonvanishing k. Indeed, below the singularity position, i.e. for k 2 < k 2 , the quantity 

1 — H§ r log 4r ) is negative. Thus, for A; 2 < fc 2 , <5(/c) becomes complex, indicating that the 



48tt 2 " 6 fc 2 ^ 

Faddeev-Popov operator (—d^D^) has ceased to be a positive defined quantity, namely one has 
left the region Co- It remains the singularity at k 2 = 0. This singularity has a simple interpre- 
tation. It means that we are approaching the boundary l\ of Co- In other words, at k 2 = we 
feel the fields on the horizon l\. 

Therefore, following pQ, a possible characterization of the factor V(Co) can be obtained by 
computing the connected two-point ghost function and by requiring that it has no poles at 
nonvanishing k. 

• Remark 

For a better understanding of the previous statement about the characterization of the 
factor V(Cq), it is useful to remind here that the region Cq is defined as the set of all 
transverse connections for which the Faddeev-Popov operator is positive definite, namely 

<rab rabc AC 



Co = {Ap, dA = 0, -dp [d^8 ab - f abc A°) > 0} . 

In the region Co, the Faddeev-Popov operator is invertible, its inverse [—0^ (d fl 5 ab — / abc ^)] 
becoming large only when approaching the horizon l\, due to the existence of a zero mode. 
Therefore, denoting by G (k; A) the color singlet Fourier transform of [—d^ (d fl 5 ab — / abc j4^)] 

xab r , \ i - l 

G ^ A ) = H -Jp^l < k I ["^ (V - f abC Al)\ | k > , (225) 

ab 

we shall require that G(k; A) has no poles for a given nonvanishing value of the momentum 
k, except for the singularity at k = 0, corresponding to the Gribov horizon l\. Finally, we 
remark that expression §225\ can be evaluated order by order in perturbation theory. It 
can be obtained by computing the connected two-point ghost function in the background of 
the gauge field A^, which plays the role of an external field. This will be task of the next 
section. 
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3.5 Characterization of V(C ) 

In order to characterize V(Cq), we start with the expression of the connected, color singlet, 
two-point ghost function 

= N j DA^5(dA)e- SYM g(x,y;A) , 

(226) 

where Q(x,y;A) stands for the connected ghost two-point function with the gauge field A a ^ 
considered as an external classical field, namely 



(227) 

We shall evaluate G(x,y;A) up to the second order in perturbation theory. Making use of the 
Wick theorem, we obtain 

(c a (x)c b (y)) c = (c a (x)c b (y) (l + J d 4 Xl (d^{ Xl )) f mnp A^)^) 
+±/d 4 a; 1 d 4 s 2 (^ 1 ^(x 1 ))/-^(x 1 )cg(x 1 )(^ 2 ^(x 2 ))/^(x 2 )c*(x 2 ))) c 

(228) 

where (%(x), c b {y) stand for free fields. Performing the Wick contraction, yields 



[c a (x)c b (y)) c = {c a (x)c b (y)) - J d^f^A^x,) (<%{x)<*( Xl )) S£ {c^(x 1 )c b (y) 

+ J d^d^r^A^r'A^) (cgOcJcgOri))^ (^( Xl )4(x 2 ))d^ (4(x 2 )c b (y) 

(229) 

From 

(<%(x)4(y)) = S ab g (x-y), (230) 
with Oq(x — y) being the free ghost propagator, we get 



c a (x)c b (y)) c = 5 ab g (x -V)-J d 4 Xl G Q (x - x^Goix, - y)f bna A^( Xl ) 
+ J d 4 Xl d 4 x 2 g (x - x 1 )d^ i g (x 1 - x 2 )d» 2 g (x 2 - y)f tna A^( Xl )f brt A:(x 2 ) . 

(231) 

Finally, for Q(x,y;A) we obtain 

g( x ,y ] A)=g (x-y)- 1 ^— [ I d i x 1 g (x-x 1 )d^g (x 1 -y)f ana A;(x 1 ) 
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(232) 

which, due to 

f aab = q ^ jacdjbcd = ^aft ^ (233) 

reads 

G(x,y;A) = Qo(x - y) 

-j^rfd^d^goix - x 1 )dH 1 g (x 1 - x 2 )d» 2 g (x 2 - y )A a l (x l )A a u (x 2 ) . 

(234) 



It remains now to take the Fourier transformation of expression (|234|) . It turns out to be useful 
to work in a finite volume V, taken here to be a four-dimensional hypercube V = L 4 , and to 
perform the thermodynamic limit, V — ► oo, at the end. 

• Remark 

The following conventions will be adopted in a finite volume. For the Fourier transforma- 
tion of the fields ip a = (A®, c a , c b ) we have 



VV q 

with 

q = — (ni, n 2 , ns, n^) , rtj integers. (236) 
lj 

The inverse of the Fourier transformation is easily obtained by making use of 

d i xe i(q-q')x = y 6 ( 23? ) 



99 
V 



Thus 



ip a (q) = [ d i x^ a {x)e- iqx . (238) 
wV Jv 

Also, 

£= E > ( 239 ) 

q ni,n2,n3,ri4 

and, in the thermodynamic limit, V — > oo, 

Let us evaluate first the Fourier transformation of the free ghost propagator 
[c? (x)c>(y)) = j DcDcl%{x)&i3,)e- * ***** 

±_J2e i{qx+py) N J DcDclf Q {q)cl{p)e-^K~ ca ^ k2ca 



(-*) 



V 

q,p 

J{qx+py) }_xabz 
y Z-^i q2° °q{~v) 

q,p 



Ir e *-j)V. (241) 

V ^ q 2 
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Therefore, for Qq{x — y) we have 

w*-v) = ^Y,j**~ v) ^- ( 242 ) 

q 

We are now ready to evaluate the Fourier transformation of Q(x,y; A). Setting 

G(k; A) = ^J d 4 xd 4 ye^ x -^g(x, y; A) , (243) 



we get 



G(k;A) = 1 I d 4 xd 4 ye lk ^G {x-y) 
V Jv 



Thus 



'v 

J d 4 xd 4 yd 4 Xl d 4 x 2 e ik ^ {G (x - x 1 )d^G (x 1 - x 2 ) 

d» 2 g (x2 - y)Al(xi)A a „(x 2 ) ) . 

(244) 



+_JL__L J2 J d 4 xd 4 yd 4 x l d 4 x 2 ( e «(*+9) e -*(*+0 e «i(p-9+«) 

qplur ^ 

x e^-^^^A^Alirfj 



N 2_l ^5 6 q-k 5 l-khp~q+u)^{l-p+r)^^ 2 j2^{u)A a v {r) 



qplur 



p 

k 2+ N 2_ 1Vk 2 k 2 2^ {p _ k) 2M P) A ^P) 



Finally 



with 



if N 1 1 ^tt - p)„ K \ 

g (1 + g (1 _ 4^ ■ (245) 



^■? H-^) ' (246) 

TV 1 1 _ (A; - g)„ fc„ 
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Summary 

For the Fourier transformation of Q(x,y; A) we have 

S{k;A) = i / d 4 xd*ye it<x -">Q(x,y;A) , (248) 

g ^ A) "h { i-o(KA)) - < 249 > 

with 

N 1 1 v^ft - q)„ K 

*(k,A) = (fc _; f W-lWto) ■ (250) 

In the thermodynamic limit, V — > oo ; 

^■^-ra/^W^- 4 ^' (251) 

Remark 

Notice also that, according to equation the expression for the Fourier transformation 

implies that the propagation of the ghosts in the field occurs with conservation of 
the ghost momentum k. 



3.6 The no-pole condition 

We can now establish the no-pole condition for the two-point ghost function. Prom expression 
(|249|) it follows that the no-pole condition at finite nonvanishing k can be stated as 



a(k,A)<l. (252) 

Moreover, following pQ, condition 1)252(1 can be simplified by recalling that in the Landau gauge 
the field A®(q) is transverse, namely 

q^{q) = . (253) 

From 

q^{- q )A a u {q) = q u A«(-q)A a w (q) = , (254) 



we can set 



AX-q)A*(q)= U (A)(6 ia ,-^) (255) 



Contracting both sides with 5n U , it follows 



u{A) = ±Al(q)Al(-q) . (256) 



Therefore, for o~(k, A) we obtain 



q 



(257) 
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In the thermodynamic limit 



a(k,A) - - n2 - ^~^2~^nu(,k) , 



= / (^(^A(-S)) (V " ^) • (258) 

As it will be checked later on, the quantity (A^(q)A1(—q)) turns out to decrease with q 2 , so 
that a(k, A) decreases as k 2 increases. Hence, as no-pole condition one can take 



a(0,A)<l, (259) 

where 

*(0.A) = 1 -J^[ . , 

(2vr) 4 q z 



This expression follows by observing that 

d A a 1 



v(o) = j ^V 2{Al{q)Al{ ~ q) A 5 " v ~ q ~f 

= S ^^ {m)A * { - q)) ' (26i) 



(2vr) 4 <7 

where the last equality follows from Lorentz covariance. In fact, setting 



and contracting both sides with 5^, we get 

j= \S^? {m)Al{ - q)) - (263) 

• Summary 

In the Landau gauge, the factor a(k, A) is given by the expression 

For i/ie no-pole condition for the two-point ghost function we have 

a(0,A)<l, (265) 

ff(M) = Ijv^i^? (jl5(9)AS( - 9)) ' (266) 

which, in the thermodynamic limit, reads 

^,A) = IJL^ J |^ (A$(q)AU- q )) . (267) 
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3.7 An expression for V(Co) 

According to [Q, the expression for the factor V(Cb) which implements the no-pole condition 
(|259|) in the path integral can be taken as 



V(C )=6(l-a(0,A)) , (268) 
where 0(x) stands for the step function^. Therefore, for the partition function Z we have 

Z = M j DA^5(dA)e~ SYM det (-9^) 0(1 - a(0, A)) . (269) 
Using the integral representation for the step function 

/ioo+e jo 

+J&'> (270) 

we arrive at the expression 

Z = M J ^DA ll 5(dA)e l3 ^°' A ^e- SY ^ det (-9„£> M ) , (271) 

which is suitable for analyzing the structure of the gauge propagator. This will be the task of 
the next section. 

3.8 The gluon propagator in the Landau gauge 

In order to work out the gluon propagator, it is sufficient to retain only the quadratic terms in 
the expression for the partition function Q271JI . i.e. we start from 



(272) 



where the limit a — > has to be taken at the end in order to recover the Landau gauge. Passing 
to momentum space, we get 

J 2mf3 

= Nj d ^ DA ^^ Aa MW-,) (273) 



with 



Thus 



Q% = ( a\ v + ( i - l) w + ^ ) 6« b . (274) 

^nuadr = / ^ (det Q%) ~ k . (275) 



^9(x) = 1 for x > 0, 6[x) = for x < 0. 
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From 

(detQ£t) 

"f(^-i)E 9 iog(« 2 +^^^), 



' 2 = g-ilogdetS^ 



. -AT/ge^), (277) 



= e V jv^-i ^ <rV ^276) 

we get 
with 

/(/?)=/?- log/3- ^(^ 2 -l)E lo gG 2 + 1^273) • ( 278 ) 
Expression 1)277(1 can be evaluated at the saddle point, namely 

^quadr ~ ^ , (279) 

where (5q is determined by the minimum condition 

/'(A)) = , (280) 

which yields 

' ~ /3b 4 F ^ 4 i ftf .9 2 J_ 
Taking the thermodynamic limit, V — > oo, and setting 



4 /3 iVff 2 1 , . 

T = iV23T2F ^->oo, (282) 



we get the following gap equation for the parameter 7 

3Ng 2 f d A q 1 



(2vr) 4 q 4 + 7 4 



(283) 



where the term 1//3q has been neglected in the thermodynamic limit, according to ea. ()282|) . The 
parameter 7 has the dimension of a mass. It is defined by the gap equation (J2HSJ)- In particular, 
from 

A 1 ^4 f A , q 3 1 . A 2 . 00/| . 

(2^7T7 = (2^ 7o dq ¥T^ = 16^ log ^ ' (284) 

where O4 = 2-/r 2 stands for the four-dimensional solid angle and A is the ultraviolet cutoff, one 
gets 

O 647T 2 1 

7 2 = A 2 e 3iv _ (285) 

To obtain the gluon propagator, we can now go back to expression (|273|) which, after substituting 
the saddle point value (3 = /3o, becomes 

iW = M 1 f DA^ VS c&W^™ , (286) 



with 

■\ab _ ( ( „2 , 1 \ x 1 ^ 1 i \ „ „ \ jfab 



QSS- = ( ( « + 72 ) V + ( - " 1 ) w» ) 5 • ( 28? ) 
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The gluon propagator is obtained by evaluating the inverse of Qz, an d taking the limit a —* 0. 
After a straightforward computation we get 

VfM = (W^Hl)) = **V^^J (<V " Qj f-) ■ (288) 
One sees that for large q, q 2 S> 7 2 , one recovers the usual perturbative behavior 

5V 1 (V - ^) • (289) 

However, for small values of q, corresponding to the infrared region, the behavior of the gluon 
propagator deeply differs from the perturbative behavior. Notice in fact that T>?^,(q) is sup- 
pressed in the infrared. 



• Remark 

In the thermodynamic limit, V —* oo ; we have neglected the term l//?o in eo. . is a 
consequence, the factor V(Co) in eq. ^221\ becomes equivalent to the 5-function 8 (1 — <r(0, A)) . 
T/ifs means that the significant range of integration in the path-integral turns out to coin- 
cide with the region near the horizon l\. 

• Summary 

As a consequence of the restriction of the domain of integration up to the first horizon l\ 
in the path-integral, the gluon propagator in the Landau gauge gets deeply modified in the 
infrared, namely 

K(q) = (A«(q)A b u (-q)) = ^V^I (<W " ^r) ■ ^ 

The parameter 7, known as the Gribov mass, is defined by the gap equation 

3Ng 2 f d 4 q 1 , , 

Wf¥T7' = u (291) 



• Remark 

Needless to say, the integral entering the equation §291\ is divergent. Within the present 
approximation, the gap equation (j£,9i|) has to be understood in a regularized way by means 
of the introduction of a cutoff A , as done in equations \28$ , Qff#5|) . In order to have a 
more precise meaning of this equation and of the Gribov parameter 7, we should be able 
to introduce a suitable set of counterterms allowing for a renormalized version of the gap 
equation (jgffij) . In other words, we should have at our disposal a local renormalizable effec- 
tive theory implementing the restriction to the region Co. Without entering in details, we 
mention that such a local formulation has been constructed by Zwanziger /?/. Remarkably, 
the resulting effective theory implementing the restriction to the first Gribov horizon l\ 
turns out to be renormalizable /?/. 



3.9 The ghost propagator in the Landau gauge 

It remains now to discuss the infrared behavior of the ghost propagator, which is obtained from 
expression \249\ upon contraction of the gauge fields, namely 

e < fc > a F(T^W) ' (292) 
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with 



a(k) - N 1 f dS {k - q) " K U*l a)A'(a)) 
a[k) ~ N*-lk*J (2tt) 4 (k-g) 2 W q)AM) 

= N 9 -75- / T-TZT; To , I ... * 



fc 2 J {2^{k-q) 2 q 4 +l 4 V " </ 2 

(293) 

Let us analyze the infrared behavior, k ~ 0, of (1 — <r(A;)). Making use of the gap equation Ifffflij) 
and of 

d 4 q 1 / q^Qv\ 3 /" d 4 <? 1 



(2vr) 4 ? 4 + 7 4 V) ~ 4 W (2^) 4 ? 4 + 7 4 ' (294) 

it follows that 

Thus, for (1 — cr(k)), we obtain 

= y (2V) 4 V ~ Jk^qf) V ~ 



(2tt) 4 V (/c — q) J <P + 7 

AT 2 /" rf4 g ( fc2 ~ 2k i) 1 ^ _ Se Qv 

" ^ 7 (27r) 4 {k _ q) 2 g 4 +7 4 \V g 
"fc2 



,2 



Atf^VW , (296) 



where 

v (k) - f A (fc2 " 2kq) 1 6 ^ f297l 

^^-y ( 2 .) 4 ( ^-,) 2 9 4 +7 4 V v " ? 2 J • (297) 

From this expression one sees that Vuvfk) is convergent and non singular at k = 0. In fact 

7V(°) = > ( 298 ) 

from which it follows that, for k ~ 0, 

?V(^o * * 2 / A^TX 1 ^ f V " ^ ■ (299) 



Since 



/ 



(2tt) 4 q 2 <7 4 + 1 A 



d*q 1 1 ^ _MA - * S /" ^ 1 1 



(2vr) 4 q 2 q 4 + 7 4 V ^ q 2 J ^ J (2vr) 4 q 2 <7 4 + 7 4 



<5 



4(27r) 4 y V + 7 4 

3 SI4 7T 1 



we get 



^8(2^) 4 2 7 2 

^iLt 2 ' (300) 



-P^^^k'S^^-^. (301) 
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Therefore 

and, for the infrared behavior of the ghost propagator 

1287T7 2 1 

^)^o-^^ F - (303) 

One sees thus that, while the gauge propagator is suppressed in the infrared, the ghost propagator 
is enhanced at k « 0, being indeed more singular than l/k 2 . 
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4 Conclusions 



• Following Gribov's suggestion, we have discussed the implementation of the restriction of 
the domain of integration in the path-integral up to the first horizon l\. As a consequence 
of this restriction, we have seen that, in the Landau gauge, the gluon and ghost propagators 
get deeply modified. The gluon propagator is suppressed in the infrared region, while the 
ghost propagator is enhanced. These remarkable features might signal that the Gribov 
copies could play a crucial role for a better understanding of the behavior of Yang-Mills 
theories in the infrared. Let us conclude this short excursion through Gribov's work by 
mentioning a few important results which have been obtained in the last two decades. 

• General properties of the Gribov region Co- 

The Gribov region Co is defined as the set of the gauge connections {A^} which are trans- 
verse, dA = 0, and for which the Faddeev-Popov operator is positive definite, —d^D^ > 0. 
The boundary of Co is the first horizon l\, where the first vanishing eigenvalue of the op- 
erator —d^D^ appears. General properties of the region Co have been established, namely 

— The region Co is convex and bounded in every direction /<§/. Essentially, this means 
that any point on l\ can be seen to have a finite distance to the origin of field space. 

— Every gauge orbit passes inside the Gribov horizon l\ JjjJ?. This result provides a well 
defined support to Gribov's proposal of restricting the domain of integration in the 
path integral to the region Cq. 

— The configuration A^ = is contained in Co . This means that the usual perturbation 
theory lies within this region. 

• Nowadays, it is known that the Gribov region Co is not free from Gribov copies, i.e. 
Gribov copies still exist inside Co ^UHl^DIS]"- To avoid the presence of these additional 
copies, a further restriction to a smaller region, known as the fundamental modular region 
A, should be implemented. However, it is difficult to give an explicit description of the 
region A. A review on the implementation of the restriction to the modular region A at 
the Hamiltonian level can be found in WI) . 

• In spite of the presence of copies inside the first horizon, Gribov's suggestion of restricting 
the domain of integration in the path-integral to the region Co captures nontrivial nonper- 
turbative aspects of Yang-Mills theories, as expressed by the infrared suppression and the 
infrared enhancement of the gluon and ghost propagators in the Landau gauge. Recently, 
it has been argued in ,13j that the additional copies existing inside Co have no influence on 
the expectation values, so that averages calculated over A or Co are expected to give the 
same result. It is worth mentioning that this behavior of the gluon and ghost propagators 
has received many confirmations from lattice simulations |141 1151 1161 1171 1181 119j . Also, 
the suppression of the gluon propagator and the enhancement of the ghost in the Landau 
gauge have been obtained within the Schwinger-Dyson approach |2U| . 

• Finally, we remark that a local action for Yang- Mills theories implementing the restriction 
of the domain of integration to the interior of the first Gribov horizon has been obtained 
by D. Zwanziger |7|. The restriction to the region Co is achieved through the introduction 
of a nonlocal horizon function in the Boltzmann weight defining the Yang-Mills measure. 
This nonlocal term may be written in local form through the introduction of suitable 

"For instance, the existence of additional copies inside the Gribov region Co can be inferred by means of 
Henyey's example as discussed in 
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additional fields. Remarkably, the resulting local action turns out to be multiplicatively 
renormalizable to all orders, obeying the renormalization group equations [7j. 
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A Appendix A. Notations 



The pure Yang- Mills action in Euclidean space-time reads 
where F£ v is the field strength 

C = d,A a u - d v Al + f ahc AlAl . (305) 

The color index a refers to the adjoint representation of a semi-simple Lie group G whose 
structure constants are given by f abc . The generators {A a } of the gauge group G are chosen to 
be anti-hermitian 

A a = -A at , (306) 



with 

Thus, 

where 



A a , X b 



f abc X c . (307) 
F„ v = d^A v - d v A^ + [Ap, A v \ , (308) 



= A«\ a , = F^ v \ a . (309) 

For the gauge transformations we have 

A^ > A^ = S^d^S + S A^S , 

S = e aaxa , (310) 



from which it follows 
At the infinitesimal level, 



F^ = S^F^S . (311) 



S = 1 + a, (312) 
a = a a X a , 

one has 

6 A,j, = d^a + [An, a] = L> M a , (313) 

SF^ = [Fn U ,a] . (314) 

In components, these transformations read 

5AI = dna a + f abc Ay = D- b a b , (315) 

5F a v = pbcpb^c (316) 

with the covariant derivative defined as 

Dn = dn + [A^ } , (317) 

D ab = 5 ab dn - r hc A C n ■ (318) 
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B Appendix B. The Gribov pendulum 

As we have seen, eq. (|4()|) . the condition for the existence of Gribov copies gives a differential 
equation corresponding to a damped pendulum under the action of several forces, see Fig. 14. 
Let us discuss here its equations of motion**. 




Figure 14: The Gribov pendulum 



From Fig. 14, we have 

d 2 y 7r 



(h , (2 + 4/ 2 )+4/ lC os(a--)+Tcos(^-a) , (319) 

dp 1 X 7T 7T 

= 4/i - Tcos(a - -) + 4/isin(a - -) , 



where T stands for the tension. From 



y = — cos(-7r — a) , (320) 
x = sin(-7r — a) , 



"It is assumed that the pendulum has unit mass m — 1 and unit length 1 = 1. 
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we have 



Recalling that 



d 2 cos(tt — a) .„ . , . . , , 7T, _ , , 
^2 " = -(2 + 4/ 2 ) + 4/ 1 cos(a--)+rcos(7r-a) 

d? sin(7r - a) . . m , 7r . . , . , vr. 
^ = 4/i-Tcos(a--) + 4/isin(a--) . 



COS(7T — QJ = — cos a , 
cos(a — —) = sma , 



it follows 



d cos a 

dT 2 

d 2 sin a 
dr 2 



sin(7r — a) = sin a , 
sm(a - — j = —cos a 



-(2 + 4/ 2 ) + 4/i sin a — T cos a 
4/i — T sin a — 4/i cos a , 



i.e. 



<i 2 cos a . . . . 2 

sma ^ — = — (2 + 4/2) sm a + 4ji sm a — i sin a cos a , 

cos a — —5 — = 4ji cos a — J sm a cos a — 4/i cos a. 

Eliminating the tension T, one gets 

d 2 cos a cf 2 sin a . 2 ,. 2 

sm a — —5 cos a — — = — = — (2 + 4/2) sm a + 4/i sm a — 4/i cos a + 4/i cos 

dr z dr A 

which becomes 

d (a! sin a) 0? (a' cos a) 

— sm a : cos a ; = — (2 + 4/2) sma — 4/i cos a + 4/i . 

ar ar 

Thus 

a"- (2 + 4/2) sin a + 4/i(l -cos a) = 0. 
Finally, adding the damping term a', one gets 

a" + a — (2 + 4/2) sin a + 4/i(l — cos a) = , 

which is precisely the Gribov condition (|40jl. 
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C Appendix C. Brief introduction to homotopy and winding 
number 

We shall give here a short introduction to the homotopy and to the winding number, following 
the references by S. Coleman j^J and by P. Goddard and P. Mansfield :22 . 

C.l Homotopy 

We shall be interested in the study of continuous mappings {(ft} between the n— dimensional 
hyper-sphere S n and the coset space Ai = G/H, where G is a Lie group and H a subgroup 

(ft(x) :S n ^M, xeS n (329) 

• Definition: Two continuous maps (ft, g are said to be homotopic if there exists a map 
F(x,t), with < t < 1 

F(x, t) : S n x [0, 1] -»■ M , (330) 
which interpolates continuously between them, namely 

F(x,0)=(ft, and F(x,l) = g . (331) 

The existence of a homotopy between (ft and g will be denoted by (ft ~ g. The set of maps 
{(ft : S n — > Ai} can be divided into disjoint classes of mutually homotopic maps, the homotopy 
classes, denoted by it n (A4). 

C.2 The winding number 

In the case in which Ai = S n , i.e. 

(ft:S n ^S n , (332) 

it can be shown that the equivalence homotopy classes are labelled by the winding number: two 
maps, (ft, g : S n — > S n , can be continuously deformed into one another if and only if <ft(x) and 
g(x) cover S n the same number of times as x covers it once. Thus ir n (S n ) = Z, the set of all 
integers. 

• Example 1. g : S 1 S 1 

Let us consider the mapping g{9) : S 1 — > S 1 , where 6 E [0, 2tt] and S 1 is the unit circle. 
For the identity map we have 

= e l(> = cosfl + isintf . (333) 

As E [0, 2tt] covers S 1 , g^(9) covers S 1 once. The map 

g (2) ^ = e i2B = cog 2Q + ■ gin 2e ^ ( 334 ) 

covers S 1 twice as 6 covers S . In general, 

g^\6) = e ive = cos u6 + ism u6 , (335) 

with v integer, covers S 1 v— times as 9 covers S 1 once. The integer v is called the winding 
number. It measures the number of times we wind around S 1 as we go once around the 
circle in two-space. 
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• The following result holds: 

Every mapping from S 1 — > S 1 is homotopic to one of the mappings g^ u \0), with v integer. 
From this result it follows that we can associate a winding number with every continuous 
mapping from S 1 — > S 1 . 

The winding number can be represented by the integral 

"=-s/„ * («"'<»>) *»""<*>■ (336) 

Moreover, it can be proven that the quantity J" 27r d8g~ 1 dgg is invariant under continuous de- 
formations, i.e. it has a topological meaning. Indeed, denoting by Sg = i(S\)g a general 
infinitesimal deformation, where S A is an infinitesimal real function on the circle S 1 , it follows 

S^T 1 = 1 => (Sg) g- 1 = -g (Sg- 1 ) , (337) 
Sg' 1 = - g ~ 1 (Sg)g- 1 = -ig^SX . 



Thus 



S r dOg^deg = r d6 (-ig-HXdeg + iff" 1 (cVA) g + ig^SXdgg) 
j o Jo 

= i / dOdeSX = . (338) 
Jo 



Also, if 



g(e)=gW(e)gM(0), (339) 
it follows that 

v = v x + ^2 • (340) 

• Example 2. 5 : S 3 -> SU(2) 

Let us discuss now the homotopy between the three-dimensional hyper-sphere S 3 and the 
group SU(2) 

g:S 3 ^ SU(2) . (341) 

The unit hyper-sphere 5 3 can be parametrized by local coordinates x^, fj> — 1, 2, 3, 4, x^x^ — 
1. We may also choose three angles {6±, 62, 63} to parametrize S 3 . The group SU(2) is 
the group of unitary unimodular two by two matrices. Any such matrix can be written as 

g = al + ib -~a , (342) 



where a, b are real parameters and I, a are the unit and the Pauli matrices, respectively. 
From 

-(!*)• °>-(-<o)- —(J -0- (343) 

it turns out 

9={ a+ h ih * h 62 + .!M , (344) 

V -fo + ih a -lbs y ' v ; 

so that the condition det g = 1 giwes 

det 5 = 1 =► a 2 + ~b ■ ~b = 1 . (345) 
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Thus SU(2) has the topology of the sphere S 3 . Therefore, the mapping Jff^iH becomes a 
mapping between two hyper-spheres S 3 

g:S 3 ^S 3 , (346) 



the homotopy classes being classified by the winding number tt^(S 3 ) = Z. Examples of the 
mapping g : S 3 —> S 3 are given by 

</°) (x) = 1 , trivial mapping v = , (347) 
g^ (x) = X4I + ix ■ ~~a , identity mapping v = 1 , 

g( u \x) = (^g^\x)J , ^integer, windings. 

• The following results hold: 

• Result 1. Every mapping from S 3 to S 3 is homotopic to one of the mapping g( u '(x) of 
equation fl^7| ). 

• Result 2 (R. Bott). Let G be a simple Lie group. Any continuous mapping from S 3 to 
G can be continuously deformed into a mapping of S 3 into an SU(2) subgroup of G. Thus, 
everything that can be established for SU(2) is true for an arbitrary simple Lie group, in 
particular for SU (iV) . 

In the case of the mapping g : S 3 — > S 3 , the expression ()336j) for the winding number generalizes 
to 

where dSn stands for the surface element of S 3 . This expression can also be rewritten as 



C.3 Application: Instantons in Euclidean Yang-Mills theories 

As an application of the homotopy and of the winding number, let us discuss here the instanton 
solution of Euclidean Yang-Mills theories. Instantons are classical solutions of the equations 
of motion of pure Euclidean Yang-Mills theories which have finite action. Let {T a } be the 
anti-hermitian generators of a Lie group G, 

T a ,T b ] = f abc T c . (350) 

Following Coleman the Cartan inner product is defined by 

T a ,T b ) = 5 ab . (351) 



For instance, in the case of SU(2), which we shall take as gauge group G, we have 

T a = - % - a a , Tr (r a T b ^j = -^5 ab , (352) 



rj-iQ, rj-ib 



1 

4 L 



a , a 



-e a = e T 



55 



Thus 

(T a , T b ) = -2Tr (r a T b ^ . (353) 
Let us start with the Yang-Mills Euclidean action 

Sym = ^ J d A x(F^,F^) , (354) 

with 

F„» {A) = d^A v - d u A„ + [Ap, A v ] . (355) 
The classical equations of motion are 

DpF^ = d^F^ + [Afj,, F^} = . (356) 

In order to have finite action, and recalling that d 4 x = dQ^dr, one requires that Fp U falls off 
faster than 1/r 3 when r —* oo, namely 

F^y ~ 1/r 3 for r oo . (357) 

This condition implies that, when r — > oo 

^ = gT 1 ^ + 0(l/r 2 ) for r -> oo . (358) 

Notice that, for a pure gauge configuration, A^ = g^d^g, one has 

F^{g- l dg) = . 

In fact 

F^g^dg) = d^dvg) - d v {g~%g) + [(jg-%g), (g~ l d v g)} 

= i 9 ^ 1 ) d u g + g^d^g - {dug' 1 ) d^g - g~ x d v d^g 

+g~ l {9^) g~ X d v g - g' 1 {d v g) g^d^g 
= {d^g' 1 ) d u g - {dug' 1 ) d^g - {d^g' 1 ) d v g + {d u g~ l ) d^g 
= (359) 

The boundary of the four-dimensional Euclidean space-time at infinity, r — > oo, is given by the 
three hyper-sphere S^. The behavior of the gauge field A^ at infinity, ea. (|358j) . allows to define 
a map between the hyper-sphere and SU{2); 

g{x) : Si -» SU{2) . (360) 

Since SU{2) has the topology of S 3 , the mapping (|36U|) can be characterized by the winding 
number v corresponding to the homotopy ^(S 3 ) = Z. This means that the classical solutions 
of the equations of motion in pure Yang-Mills with finite action can be classified by the winding 
number v. 

In order to find classical solutions of the equations of motion, it is useful to consider the identity 

J d A x {F^v =F* F^F^ =f* F^) > , (361) 



where *F tlu is the dual of F^ u 



Fpv — ^^lupaFpa , (362) 

£^vpu£^uX5 = 2 {5 p \5 a s — 5 p s5 a \) . 
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From eq. ljUfilj) we have 

J d 4 x (F^ =F* F^F^ =F* Fp,,) = J d 4 x ((F^,F^) + (*F^,* F^) T 2 (F^* F^)) (363) 



Since * = -F^, it follows 



The bound (|364() is saturated when 



-L / d A x{F^,F^) > -f- 



4<? 2 



/" <i x (i^y, F^y t 



(364) 



F^jjy — i i^jjy — ^^pvpa^pa • 



(365) 



This condition is a first order differential equation. The solutions to the self-dual equation 



1 



£ ' ti.unn-Fn, 



(366) 



are called instantons (anti-instantons are solutions of Fn V 
we have the equality 



t F flll ). Thus, for an instanton, 



pv ; Fpv ) 



J d x (F^ U: F^v) 



(367) 



It is useful to observe that, from the self-dual condition, F„ u =* Fn U , one has 



DpF^ v — r.^pvpa-D^Fpa — , 



(368) 



due to the Bianchi identity. This means that instantons are solutions of the equations of 
motion of pure Euclidean Yang-Mills theories. 

Another important property of the instantons is that they give vanishing contribution to 
the energy-momentum tensor Qp U , as it is apparent from 



®^ ~~ 2g 2 + * Fpp * * F " p ^ 



(369) 



It is useful to show now that the quantity f d A x (F^* F^) = f d A x (FjZ v *F®) is directly related 
to the winding number v. In order to establish the relationship between f d 4 x (Fp U * Fp V ) and 
v, let us first prove the identity 



Tr (F^*F^) 

Qp 



dpGp , 



(370) 



Tre 
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Making use of the Bianchi identity, Ep U p a DpF pa = 0, we have 

d/iGfj, = Tre^d^vFp,, - 2 -A v A p A^ (371) 
= Trs^pu {{dpAy) Fp a + A u dpFp a -2(dpA u )ApA <T ) 
= Trs^pu ^— {dpA u — d u Ap) Fp a — A v [Ap, Fp a ] — {dpA u — d v Ap) ApA^ 

= Tr£p U p a ^—Fp U Fp a — — [Ap, A u ] Fp a — A v ApFp a + A V F p a Ap Fp U ApA a ^j 

-\-Tte p^pa {[Ap,, A u ] ApA a ) 
= Tr£p U p a ^—FpjjFpu ApA u Fp a + ApA u Fp a + ApA v Fp a ApA a Fpy^ 

^FrEpypu i\Ap, A u \ ApA a ) 
= TrEp U p a ^—Fp^Fpu + [Ap, A v \ ApA,j^j . 

The term TrEp Upa [Ap, A v \ A p A a vanishes due to the Jacoby identity 
FrEpypu [Ap, A u \ ApA a = —TrEpjypa [Ap, A v \ [A p , A a ] 

a a Ab ac Ad fabm fcdnrp rpmrpn 

— 2 i 1 "^ V v P ■> 

1 Ad *b ac /id fabm fcdm 

— -^fJ-upcr^p^u^-p^al I 

= . (372) 

Thus, we conclude that 

dpQp = TTEp V p a l -Fp V F pa = Tr {Fp V *Fp U ) . (373) 
Moreover, making use of the Stokes theorem we obtain 

Tr f d 4 x {Fp U *Fp U ) = J d 4 xdpGp = Tr dSpEp Upa (a u F P(T - ^A v A p A^j . (374) 

For a classical solution of the equations of motion with finite action, we have that on the hyper- 
sphere at infinity 



Ap = g-'dpg on Si (375) 
Fp 

Therefore, for an instanton solution 



Fpu = on . 



Tr J d A x{Fp V *Fp V ) = Tr dSpEp Vpa ^A v F pa -^A v A p A a 

= ~~ —Tr j dSpEp V p a {ApApAfj) 
<3 Js' A 

• / dSpEp Vpa {g~ l d u g) {g^dpg) {g~ l d a g) 



(376) 



= — Tr 
3 



= -16vr 2 i/ , 
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namely 

v = -j^Tr J d 4 x (F^F^) , (377) 

or 

U= 3^j d4x {F »"* F ^ = 6^ / *™i>»t"'K»K° ■ ( 378 ) 

The expression for v, ea. (|378|) . is called the Pontryagin index, whilst the integrand is the Pon- 
tryagin density. Thus, for an instanton with winding number v 

Sym = 4^2 / ^xF^ = ^ M . (379) 

In the case of SU(2), the explicit solution for the instanton with v = 1 has been given by Belavin, 
Polyakov, Schwartz, Tyupkin, and reads 



A, = f(r 2 )(gV) 8/), (380) 



(!) x 4 + i x ■ a _ 

r 2 

f(r 2 ) - 



r 2 + p 2 



where p is an arbitrary constant, called the size of the instanton. 
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D Appendix D. Polar coordinates 



Let us remind here some useful relationships in polar coordinates: 



x = r sin 9 cos <p , 
y = r sin 9 sin <p , 

z = rcosO. (381) 

Thus, for the orthonormal basis (~e r , ~e$, ~~e ip) we have 

e r = sin 9 cos p e x + sin 9 sin p>~e y + cos 9 e z , (382) 
~~e q = cos 9 cos p~e x + cos 9 sin p~e y — sin 9~e z , 
~e v = ~e r x ~eg = — sin p~e x + cos 99 "e^ . 

Let ~v be a vector. We have 

1/ = v x ~e x + fyl^y + ^~ez = v r ~e r + ue'e'e + v (p ~e v , (383) 



with 



u r = ~v ■ ~~e r = v x sin 9 cos 99 + v y sin sin p + v z cos , 
vg = 1? ■ ~e ' q = v x cos 9 cos 99 + f y cos sin ip — v z sin , 
v v = v ■ e v = — v x s\nip + v y cos 99 . (384) 



Also, 



V J = a~ e r + e H ^« _e V' 

ar r at/ r sin 



72* 1 9 f 2df\ , 1 d (_,__ n df\ , 1 a 2 / 



V 2 / = +TT73ih^ + 



r 2 <9r \ <9r / r 2 sin 9 89 \ 89 J r 2 sin 2 dip 2 ' 
= 1 d(r 2 v r ) | 1 9(smM | 1 dv^ 
r 2 dr rsin9 39 rsin# dp 
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